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Abstract. The ground state and the excitation spectrum of strongly correlated electrons in 
quantum dots are investigated. An analytical solution is constructed by exact diagonalization 
of the Hamiltonian in terms of the A''-particle eigenmodes. 



1. Introduction 

Spatially confined charged particles are attracting increasing interest. In particular, quantum 
confined semiconductor structures allow to design materials with novel properties, see e.g. Ref. 
^ for an overview. Of special interest are quantum dots which allow to confine charged particles 
in all three space dimensions, e.g. [2]. Besides their application in nanoscale transport and 
optics, quantum dots allow to create strongly correlated many-electron states including Fermi 
liquid and Wigner crystal behaviour [3] which are of fundamental interest in many fields. While 
the ground state of these systems can be studied theoretically from first principles, e.g. using 
quantum Monte Carlo methods, these methods fail to deliver the energy spectrum. Therefore, 
various analytical concepts have been considered which allow to treat small systems of confined 
electrons. The problem of the energy spectrum of two electrons in a harmonic trap has been 
solved in Refs. [HIEI. 

Here, we extend the analytical approach of Ref. 7 to the general case of strongly correlated 
particles (e.g. electrons). We demonstrate in Sec. 1 that, in the strong-coupling limit, the A^- 
particle problem can be reduced to a superposition of (A^ • d) uncoupled harmonic oscillators. 
These results are illustrated for two electrons in Id and 2d in Sec. 3 and 4 where we demonstrate 
how to analytically construct the ground state wave function and energy spectrum in the whole 
range from strong to moderate coupling. 

2. Energy spectrum of N strongly correlated electrons 

The Hamiltonian of Coulomb-interacting electrons of effective mass m trapped in a d- 
dimensional harmonic potential is given by 

N -Ji ( N N \ 

i^ = r(p^...,p--,) + m,...,r^) = E£^+ E?-'^t+E^ ' (1) 

i=\ \ i=l i<j j 

2 

where u; is the confinement frequency, and a = ^^^^^^ is the Coulomb interaction strength, with 
eb being the background dielectric constant. 



The degree of electron-electron correlations is characterized by the dimensionless coupling 
parameter Q jSj, which is the ratio of the characteristic Coulomb energy Eb to the confinement 
energy hw. 

where Eb = = ^f^i = \J ^ is the characteristic confinement length, and qb denotes the 
effective Bohr radius. In the strong-coupling limit, Q ^ 1, the electrons are spatially localized 
around their classical equilibrium positions fo = {fio, ■ ■ ■ ,rjvo), and the confinement potential 
\^(r) can be well approximated harmonically around fo. 

V{f) = V{fo) + ^(r - fofA (r - fo) + 5V , (3) 

where Aij = ^q^q^^ with f= (^i, . . . , rTv-d) is the • d-dimensional Hesse matrix, and the 
anharmonic correction 5V is neglected, it will be restored in Sec. 5. Further, since matrix A 
is symmetric there exist square • d-dimensional matrices B and U , such that A = U'^ B U, 
where B is diagonal and U is unitary. Introducing ^ into yields 

H-V{fo) = ^{Uf)f{Uf))+l{U{f-fo)fB{U{f-fo)) (4) 

= -L{Pf(P) + hRfB{R), 
2m 2 

where p = {pi, . . . ,PN-d), and the new momenta and space coordinates, P and R, are defined by 

P = Up, R = U{f-fo). (5) 
With the diagonal elements Ba = milf , Eq. Q simplifies to 

^■(^ / p2 \ 



i=l 



and the A'^-particle problem has been reduced to a superposition of {N ■ d) uncoupled harmonic 
oscillator modes with the eigenfrequencies 




f),(u;) = W^. (7) 



This description in terms of normal modes arises in the harmonic approximation of the 
potential around the classical equilibrium position ro, which is exact for strong electron-electron 
correlations, Q — > oo. The eigenfunctions of the A^-particle Hamiltonian of Eq. ^fi{R), can 
be written as a product of (A^ • d) one-dimensional oscillator eigenfunctions 

N-d / 52 \ 

^fi{R) =\{^nM) > i^nARi) = Cn. ^Xp --^ Hn^Ri) , (8) 



where n = (ni, . . . , nAr.rf), Ui = 0,1,2,... , and HmiRi) are the Hermite polynomials. The 
corresponding excitation spectrum around the classical ground state F(ro) is given by 



N-d 

(10) 



i=l ^ ' 



In the following two Sections, our general approach is illustrated for N = 2 strongly correlated 
electrons in Id and 2d. 

3. Two electrons in a \d harmonic trap 

For two electrons in a one-dimensional harmonic confinement 

^ = m,p.) + l^(x„..) = ^|^+(^^.\^ + ^^^| , (11) 

«=1 \ i=l i<j 

V{x)^V{xo) + ^^ {x-xq)'^A{x-xo) , (12) 

where x = {xi , X2), and xq = (^^lOj 2:20) is the classical equilibrium position of the electons, which 
follows from V V{xq) = : 

1 f 2a 

The equilibrium potential energy is 

V{x,) = l[^aWf" = nu:^^Q, (14) 
and, using H13() . the matrix A becomes 

^ V ' 

=B 

The unitary matrix [/ (jlSp defines the canonical transformation to the new variables (xi,pi) 
and (x2,P2), 

V ^2 - ^ y ^2 V ^2 - a;i - xo y ' ^ ^ 

\ =u(P' \ =^(P'+P' \ . (17) 

P2 y V ?'2 y ^2 V ^"2 - y 

In the harmonic approximation {6V = 0), the electron dynamics then can be described as a 
superposition of • d = 2 uncoupled harmonic oscillator modes 

2 / ~2 \ 

-F(xo) = 5^(1^ + ^1^2 ~2\ ^ Oi=w, Q2 = ^/3u;, (18) 

The first mode of frequency lv corresponds to the oscillation of the electron center-of-mass, 
whereas the second (breathing) mode, see Ref. of frequency \/3a; describes anti-phase 
oscillations of the two electrons. Replacing pi in Eq. (|18|) by the momentum operator —iKVxi 
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Figure 1. Ground state correlation function 
Kofl{xi, X2) for coupling parameter Q = 50 
(oowTOOas). In xi- and X2-direction one 
recognizes the different widths of the peak, 

5(a). 



and introducing dimensionless coordinates xi and X2 by (cf. Eq. Q) = oq Xj , yields the 
Hamiltonian 

H - V{xo) 



The eigenfunctions of Eq. (jT!?]) . ^'^(xi, j;2)5 are of the form of Eq. with fJ = (ni,n2), 
rij = 0, 1,2, . . . , i.e. ^'^(xi,X2) = -0^1(^1) • tpn2{x2), where 

Ipn^iXi) = Cm exp (^-^^^^i?n,(^/^Xi) , Cn, = ^^^1— y (^) * • (20) 

The two-particle excitation spectrum is then given by 

En„n, = V{xo) + e« + e^L'^ , e« = hQ, + . (21) 

The properties of the ground state wave function, ^(^q q^{xi,X2), can be illustrated using the 
two-particle correlation function in real-space 



K^{x,x')= j dxi J dx2\'^nun2ixi,X2)\'^ ■ [pix)p{x') - 5{x - x')p{x)] , (22) 
where p{x^'^) = Yll^i 5{x^'^ — Xi). Eq. (|22|) simplifies to 

Kri{xi,X2)=2\iJn^{xi)\^ |V'n2(^2)|' , (23) 

where from Eq. (ITBl) it is xi = ^7±^ and X2 = ^"^-^-^o ^ ^pj^g density plot of Kqq{xi,X2) is 
presented in Fig. 1. The ground state correlation function shows a Gaussian-like peak around the 
position {xi,X2) = (— ^,^), which has two different characteristic widths (5(0j) in directions 
xi and X2 of the corresponding eigenmodes. 

Next, consider the ground state electron density /0(o,o)(^) defined by [H] 

Pn{x) = / dxi / dx2\'^n{xi,X2)\'^p{x) , (24) 
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Figure 2. Ground state electron density 
Pofl{x) for Q = 25 (dashed curve) and Q = 
125 (solid curve). 



where for ni = n2 = we obtain 



Po,o(.x) = — { e 
Po 



+ e 



2{i + V^)h 



"v/S m uj X 



(25) 



with K being the dimensionless variance. The factor 1/po ensures that the electron density is 
normalized to the particle number, J pQfi{x)dx = 2. po.o is shown in Fig. 2 for parameters Q 
in the moderate and strong coupling limit. The variance k determines the width of the peaks 
around the equilibrium positions — ^ and Using ((T^ we obtain k oc cj^/^ oc ^. In the 

strong-coupling limit, it is k ^ 1 corresponding to strong electron localization. For Q<1, 
strong electron overlap is observed, which, however, is beyond the harmonic approximation of 
Eq. (inj, and one has to account for anharmonic corrections 5V, see Sec. 5. 



4. Extension to two electrons in 2d 

Here, for the case of two harmonically confined electrons in 2d we analyze the correlation function 
K{f,f'). Using the notation r = (ri,r2) = ixi,yi, X2,y2), the Hamiltonian reads 



1 

2m ' 



a 



In - 



(26) 



where p = {pi,P2)- The classical equilibrium distance of the electrons is given by 

= 23ao ^/Q . 



/ a 



(27) 



In the strong-coupling limit, the total potential can be expanded as in Eq. (j^J with 6V = 
resulting in 
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V2 
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B = muj'^ 
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(28) 



With unitary transformations of Eq. © the two-particle problem, Eq. H26() . is reduced to 
four uncoupled harmonic oscillators corresponding to two centre-of-mass modes ^2,3 = "^i one 
breathing mode Qi = \/^uj and the trivial rotation of the whole system with ^4 = 0. 



The ground state correlation function ((22} is evaluated to 

K,{r,r') = ^ exp (^-{Dl + Dl)-^[r- 2\ (29) 

where Dx = ^2a^' ~ ^^o^' ^ ~ (^l^i ~ ^^2! and kQ denotes a normalization constant. 

In Fig. 3, Kq is shown for the strong coupling case Q = 250 and one electron being fixed in 
three different positions. Fig. 3 c) shows the correlation function when one electron is situated 
in its classical equilibrium position, e.g. no = (^^)0). The density of the second electron now 
shows a Gaussian-like peak around its respective equilibrium position r20 = —Ho- As in Sec. 3 
(Fig. 1) one recognizes characteristic widths in different directions corresponding to the different 
eigenmodes. 



5. Anharmonic corrections 

So far, we confined ourselves to the strong-coupling limit, Q » 1, where 5V = is adequate. 
In contrast, at weak coupling, in Eq. (jSJ the anharmonic corrections 5V have to be taken into 
account. Since the derived eigenfunctions ^fi{R) =: \n) (jSJ of the harmonic approximation form 
a complete basis, the higher order terms 6V can be expanded in terms of these eigenstates. 
Here, we limit ourselves to spinless particles in order to demonstrate the general scheme. In 
the following, the ground state and the excitation spectrum for 6V 7^ is derived for the one- 
dimensional case of Sec. 3. Using the transformations (jlbj) . we introduce new coordinates xi^2 
given by xi^2 = + "^(^1 "F ^2)5 in terms of which the total potential energy has the form 

V{i) = V{X,) + Vix2) = ^Xl + ^(Xo + iX2f + , " _ , , (30) 

where x = {xi,X2) and = ^/2aQ. Obviously, anharmonic effects arise exclusively from V{x2)- 
The expansion of V{x2) around X2 = yields 

V{X2) = v{xo) + ^h^xl + ^^^$r^ ■ (31) 

j>3 ^0 
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Figure 3. Ground state correlation function KQ[r,f') for Q = 250. One electron (dot) is fixed 
in different positions n = (xi,yi) and the density of the other electron is shown, ■61. Lengths 
are displayed in units of oq. Configurations of single figures: a) ri = (0,0), b) fi = (0.8,0) and 

c) ri = ( 2^'0) — (9.96,0), corresponding to the classical equilibrium position. 



Expanding 5V in terms of Hermite polynomials Hk{x2), the sum over i in Eq. (|31|) can be 
rewritten as 

6V = ^iyk{Q) Hk{x2) , (32) 

fc>0 

where the expansion coefficients i^kiQ) are easily obtained by comparison of coefficients, and the 
Q-dependence arises from the Q-dependence of xq, see Eq. (|13j) . The Hamiltonian corresponding 
to coordinate X2 is then given by 

-J— = "2 + + L ^ ^^=(^2) , (33) 

fc>0 

where we write H = H{xi) + H{x2) + V{xq). To numerically determine the ground state 
energy and relative wave function, ^q{x2), i.e. the eigenfunction of H{x2) corresponding to 
H{x2) ^q{x2) = s <l>o(2;2), we limit the summation in Eq. to finite fcmax > 3. Then the 

ansatz reads 

m 

^o{x2) = Y,Cn{Q)-\n) , (34) 

n=0 

where \n) = ipnix2) are the basis eigenstates, see Eq. (j'iOp . Using the explicit form of H{x2), 
the coefficients follow from multiplication with (n'| : 



n=0 

.(2) 



J] CniQ) [e^^^Sn'n + (n>y|n)j = e Cn'iQ) , (35) 

n=0 

where en' is given by Eq. 1)21(1 . The matrix elements are computed as 

^max ^max 

{n'\5V\n) = ^fe(Q) (n'|Ffc(x2)|n) = ^ z>fc(Q) (n'|/7,,(y)|n) . (36) 



k>0 k>0 



where ^X2 = y, V = ^ = V^) and |n') = ipn'{x2/ ^/v)- The integral (n'|i7fc(t/)|n) is evaluated 
exactly with s = ^{n' + k + n) and Cn given by Eq. (|20|) : 

{h'\Hk{y)\n) = ^/r] Cn'Cn j dy exp{-y^) Hn'iy) Hkiy) Hniy) 

= \A/ C„'C„ < 



2''n'!fc!n!0F 



(s-n')! (s-A:)! (s-n)! 

, otherwise 

Below, we present results for the ground state wave function and the excitation spectrum 
for fcmax = 6 and different coupling parameters Q <1. The matrix eigenvalue equation H35() 
has been solved numerically with variable upper summation limit m yielding a fixed precision 
of the energy eigenvalues, Ae ~ W^^ hu>. The expression for 5V in Eq. ()31|) in terms of y 
and Q, and the coefficients i^kiQ) of Eq. ()36j) are given in the Appendix. In the following 
figures, the results for different values of the coupling parameter are shown. As one sees in the 
coefficient spectra of Fig. 4 to 6, mainly the even states contribute in the expansion of Eq. 
()34|) . Particularly, for Q < 0.1, even coefficients are several orders of magnitude larger than 
the odd ones. Also, the spectrum gets broader for decreased coupling, and the values drop off 
approximately exponentially with n. ! This property is marked by the straight lines in Fig. 4-6, 
corresponding to exponential fits according to 
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Figure 4. Square of the expansion 
coefficients Cn{Q) for Q = 1. Tlie ground 
state energy eigenvalue is e = 1.1539 huj. 
The coefficients of Eq. (50) are /(^) = 1.87, 
/(°) = 1.30. 
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Figure 5. Same as Fig. 4 for Q = 0.1. Figure 6. Same as Fig. 4 for Q = 0.01. 
£ = 2.7703 hu. f^""^ = 0.66, /(°) = 0.33. e = 8.4996 hw. /(^) = 0.2, = 0.06. 



n{e,o) 



oc exp 



(37) 



where (e, o) stands for even and odd states, respectively. The values of f^^^°\Q) are given in 
the captions. In Fig. 7, the Q-dependence of the ground state and the first and second excited 
state is shown. 

In Fig. 8, the fit parameter /^^'"^Q) in Eq. ()37|) is shown as a function of the coupling 
parameter. The Q-dependence is approximately given by: /(^) oc ^/Q and /(°) oc Q^. Thus, for 
Q < 10 the ground state wave function can be well approximated by 



^o(xi,X2) 
^o(Si) 



^0(2:^2) = ^ exp(-Ai \/Qn) ■ \n) + ^ exp(-A2 



nj ■ \n) 



(38) 



n odd 



where ^o{xi) = ^po{xl) as in Eq. ((^ . and Ai, A2 are given in the caption of Fig. 8. 



6. Discussion and Outlook 

In this paper, we presented analytical results for the wave function and energy spectrum of 
strongly correlated electrons. The results are exact for harmonic interaction as well as for 
arbitrary interaction in the strong coupling limit, Q — > 00. 

For quantum dots based on a GaAs/ AlGaAs heterostructure (e^ = 12.9, m = 0.067 me) with a 
characteristic intrinsic confinement energy of /iw = 5 meV, we obtain the confinement frequency 
u! = 7.6 THz corresponding to a coupling parameter Q = 1.3, an equilibrium distance (|13j) 
xq = 21.7 nm, and a width of the Gaussian peaks in the ground state electron density of Fig. 2 
of "v/k = 1.3xq. a situation of strong coupling of e.g. Q = 50 would require a shallow confinement 
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Figure 7. Influence of the anharmonic effects 
on the excitation spectrum: Q-dependence 
of the three lowest energy eigenvalues corre- 
sponding to H{x2). The straight lines show 
the equidistant spectrum en of Eq. (|^T|). 
valid for Q — > oo. 
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Figure 8. Q-dependence of the parameters 
j^q_ (37). The errorbars result 
from curve fitting with Eq. (37) and mark 
lines corresponding to minimum and maximum 
gradient, respectively. The lines are best fits 
with Ai = 1.81, A2 = 1.05. 



potential of huj = 0.088 ^eV leading to xq = 32.1 /im and -^/k = 0.2 xq. Stronger coupling can 
be achieved in materials with strong Coulomb interaction. E.g. for a ZnSe structure (e^ = 9.1, 
m = 0.21me) we obtain Q = 2.4 (at huj = 5 meV). 

For finite coupling, Q > 1 , an analytical result is obtained by expanding the anharmonic 
corrections of the potential energy in terms of the A^-particle eigenfunctions of the strong 
coupling (harmonic) limit. For moderate and small coupling, individual electrons overlap 
significantly. Therefore, in this case the resulting A^-particle wave function has to be anti- 
symmetrized. The corresponding generalization of our method will be presented elsewhere. 
Furthermore, there are two aspects of our approach that will be examined in future work: 

1. First principle spectrum from Quantum Monte Carlo data. Path integral Monte Carlo 
(PIMC) simulations allow to compute the density matrix of the A^-particle problem from first 
principles [31 El- However, they cannot directly yield the energy spectrum and corresponding 
wave functions. In the limit of strong correlations, the idea is to reconstruct the energy spectrum 
and wave functions from PIMC ground state results for the density distribution or correlation 
function using an expansion in terms of the derived oscillator eigenmodes 

i=l \n=0 / 1=1 L/n=0 ^" (Vj 

E.g. a fit to the PIMC the density distribution 

PQMc{r) = I d\i . . . d\N \'^o{R)? P(r) , (40) 

with R = R{fi, . . . , rXr) and p{r) = X^i^i ^{r — ri) then yields the expansion coefficients. 

2. Application to Optics and Quantum Transport. In the strong coupling limit, the A^-particle 
Hamiltonian ^ including external fields, SHext, can be written in second quantization as 

N-d . . 

H - Vifo) = J2^^i(4ai + ^j + SV + 6Hext , (41) 



where cij are expected to be bosonic field operators of tfie derived eigenmodes satisfying tfie 
commutator relation 

■ ' = <5mn ■ (42) 



In terms of these field operators, one can define Nonequilibrium Green's functions by 

9^nM = j^{&iit')^m{t)) , glnM = j^{&m{t) aUt')) , (43) 

and introduce density matrices 

Fmn{t)=ihg<^{t,t) . (44) 

The subject of ongoing work is to derive equations of motion for g^^{t,t') and g^^{t,t') or 
Fninit), respectively, the solutions of which give direct access to the optical and transport 
properties of the strongly correlated confined A'^-particle systems. 

Appendix. Explicit structure of 6V in Sec. 5 

SV in Eq. (37) can be rewritten in terms of y = ^/fjx2 and Q as follows {imax = 6) 




26 \ 26 o 23 , ^ . 23 

VvQ J 34VQ 34 g 2 3V3Q2 



Expanding 6V in (A.l) in terms of Hermite polynomials Hj,{y) yields 

k=0 

where the coefficients i^kiQ) are given by 

_ 5 • 2I V3 + 6 • 2I Q , 5 • 33 V2 + 6 • 2§ 3i g 

"°^^^ = 24g^ ' = ^1^1 ' 

, 5-2iV3 + 4-2ig , 5^23! +3-2§3ig 

'^^(^) = leg^ ' "'^"^^^ ' 

, 5-2i V3 + 2-2^g , 1 - 1 

96g2 48^234 g2 9623 V3g2 
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